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Abstract 

The superfluid phase transition of the general vortex gas, in which the cir- 
culations may be any non-zero integer, is studied. When the net circulation 
of the system is not zero the absence of a superfluid phase is shown. When 
the net circulation of the vortices vanishes, the presence of off-diagonal long 
range order is demonstrated and the existence of an order parameter is pro- 
posed. The transition temperature for the general vortex gas is shown to 
be the Kosterlitz — Thouless temperature. An upper bound for the average 
vortex number density is established for the general vortex gas and an exact 
expression is derived for the Kosterlitz — Thouless ensemble. 



1. Introduction. 

In 1973 Kosterlitz and Thouless published the first of two papers [1], [2] 
describing a phase transition in two dimensions. Four years later Nelson and 
Kosterlitz [3] used this theory to predict the celebrated universal jump in the 
superfluid density of '^He films which was subsequently verified by Bishop 
and Reppy [4] in their famous torsion pendulum experiment. The corner- 
stone of this theory is the presence of a neutral "gas" of point vortices in the 
film which interact through a coulombic hamiltonian. A natural analogy is 
made with a two dimensional gas of point charges and the phase transition 
is thought to be due to the binding and unbinding of dipole pairs. The 
high temperature phase consists of free vortices each having a circulation 
or "charge" |g| = 1 while the low temperature phase consists of oppositely 
charged vortices bound in dipole pairs. It is thought that the presence of 
free vortices breaks any long range order in the system [5] and, from the 
Hohenberg-Mermin- Wagner (HMW) Theorem [6] , [7] , we know that the ^He 
film will behave as a normal fluid. When the temperature is lowered below 
the transition temperature the gas of vortices condenses into dipole pairs 
and are no longer free. Long range order is thereby possible. Renormaliza- 
tion group techniques are then used to analyze the critical properties of the 
system and give an operational definition of the superfiuid density. 

We shall analyze the phase transition in a different, somewhat more 
traditional way. Unlike Kosterlitz and Thouless we shall not begin by differ- 
entiating between "free" and "bound" vortices but will instead treat them 
all on an equal footing. We shall apply an infinitesimally small external fiow 
or "electric" field Ej to the system which may be slowly varying in time. It 
will serve two purposes. First, it will polarize the system and, as it will be 
coupled linearly to the polarization vector Pj = Yli^otfJ^ the grand partition 
function becomes a generating functional for correlation functions. Second, 
noting that Pj is proportional to the total superfiuid momentum density g^j^ 
in the film, by looking at the superfiuid momentum-momentum correlation 
functions (gj^gf/) we can make contact with the standard two fiuid model 
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in three dimensions. As usual, the normal density will be proportional to 
the trace of this matrix and will be required to be finite in the static limit 
Ej 0. The signature of a superfluid state will be when the off-diagonal 
terms of the matrix becomes infinite in this limit, signifying off-diagonal 
long range order. The coefficient of the off-diagonal terms will then serve as 
our order parameter. 

Our approach bares some resemblence to Linear Response Theory. Tra- 
ditional Linear Response Theory as it is applied to a classical ensemble, 
however, requires a perturbative expansion in the time derivative of the po- 
larization vector Pj about a small applied field and then letting that field 
go to zero. For the vortex gas the polarization vector is a constant of the 
motion and Pj vanishes identically. The perturbation procedure ends be- 
fore it begins. Fortunately, enough information can be obtained from the 
grand partition function to demonstrate the presence of long range order as 
well as the existence of an order parameter without resorting to Perturba- 
tion Theory. Moreover, we find that we need not require that the vortex 
charges be only ±1 and will instead begin by considering the more general 
system in which the charges may be any non-zero integer. Nor shall we 
require that the net charge Q of the system to be zero, although we will 
consider the cases Q ^ and Q = separately. Our fundamental premise is 
that a superfluid state does exist in two dimensional ^He films, and that the 
Kosterlitz — Thouless — Nelson theory is the correct description of it. Our 
results will, however, require a re-interpretation of the transition itself. 

A few words about terminology and notation. By the General Kosterlitz 
— Thouless ensemble, or simply the general ensemble, we mean the ensemble 
in which the vortex circulation may take on any non-zero integer value. The 
thermodynamic average of an object O which is calculated using the grand 
partition function for this ensemble will be denoted by (O) and is a 

function of the applied field Ej. What will be of particular importance is 
the value it takes in the static limit Ej — > which we will denote by the 
subscript zero: (C)o- We may also calculate the thermodynamic average of 
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O when Ej = identically and this will be denoted by {O)\Ej=0- We shall 
find that for some O, {O)\Ej=0 7^ (C')o- When we restrict our attention to 
the Kosterlitz — Thouless ensemble, in which all the vortices must have unit 
charge and the net charge of the system must vanish, we shall denote all 
thermodynamic quantities calculated using their grand partition function 
with the superscript kt. Greek indices will run from 1 to N, the number of 
vortices in the system, while latin indices will run from 1 to 2, the dimension 
of the space. The summation convention will be used throughout this paper 
when latin indices are repeated. It will not be used for greek indices. 

The rest of this paper is organized in the following manner. Section 2 
provides a brief review of vortex dynamics in two dimensions and gives the 
motivation for our construction of the grand partition function. In Section 
3 we shall consider the case in which the net charge does not vanish and 
show that for this system there is no superfluid phase. The vortices will 
always behave as a normal fluid. In Section 4 we shall consider the case 
when Q = and demonstrate the existence of a superfluid phase transition 
at the Kosterlitz — Thouless temperature Tkt- An order parameter for the 
system is proposed and its' relationship with the average density of vortices 
is derived. All of this analysis is done away from the transition temperature 
and it is not until Section 5 that the behavior of the system at Tkt will be 
addressed. Concluding remarks may be found in Section 6. 

2. Background. 

We begin with a brief description of point vortices in Bose liquids. Let ip 
be the microscopic complex scalar field for the ^He atoms. The microscopic 
current density is 

where m is the mass of the "^He atom. We define the superfluid velocity as 
v^^ = p~^j[8] where p = IV'P is the two dimensional density of the "^He film. 
v^f is, as usual, proportional to the gradient of the phase of ip. As such, for 
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any closed path 7 in the film, 

jj-f . dl = q— (2) 
m 

where h = 2'Kh. Because i/j is a Bose field, q may be any integer. When q is 
non-zero, we say that there is a vortex with circulation q somewhere within 
7- 

At times we shall find it more convenient to work in complex coor- 
dinates. We define z = {xi + 1x2) /V^ with the corresponding definition 
v^f = {v^-^' — iv2^)/\/2. As usual, complex conjugates are denoted by a bar. 
If we now consider a vortex with circulation q^. located at a position z"-^ in 
the film, (2) becomes the contour integral 

/ v^Jdz = q^- (3) 
J-y Tna 

*^ let 



around any closed path 7q, in the complex plane encircling z°^ . As such, the 
superfiuid velocity for this single vortex is 

<^ = W^- ^ + ^« (4) 

where is any holomorphic, or meromorphic function of z with poles of 
order greater than one. Since we are primarily interested in the first order 
pole of v^f, we shall set v'^ = 0. We caution the reader that our definition 
of the superfiuid velocity differs somewhat from that used by Nelson and 
Kosterlitz [3] but agrees with Minnhagen and Warren [8] . Letting po be the 
spatial average of p, we define a total superfiuid current density j^-^ for N 
vortices located at positions {z"} with circulations {qa} by 

, N 

27rim ^ z - z^" ^ ' 

OL = \ 

as well as a superfiuid "momentum" density g^^ = mj^^ . 

Looking back at the original proof of the HMW Theorem, we see that 
in addition to Bogoluibov's inequality the conservation equation 

(a) 



corresponding to the continuous symmetry (in our case a U(l) gauge sym- 
metry) was also used. Because — log\z — w\/27r is the Green's function for 
the Laplace Operator in two dimensions, we find that ^ • J*-^ = 0. As such, 
the HMW Theorem does not take into account vortex excitations in the 
fluid. 

Since ^•J*'^ = 0, while still satisfying (2), the vortices behave as though 
they are point vortices in an ideal liquid. Treating the vortices as though 
they were point particles in and of themselves, their lagrangian is: 

N 

dt 



= -iY^q^z-'^ - IC (7) 



where 



and 



The reader is referred to either [9] for a standard hydrodynamical derivation 
of (7) or [10] for a derivation from the microscopic '^He lagrangian. e^t is 
the energy scale for the system which we have taken to be poT^Ti^ / (2m). We 
caution the reader that our potential U differs from that used by Kosterlitz, 
Thouless and Nelson by a factor of 2. As usual, a hard core cutoff has 
been introduced at a distance a to prevent any infrared divergences in the 
partition function. 

Taking our generalized coordinate to be its' canonical momenta is 
then —iqaz" and we see immediately that £v is linear in the momenta. 
Consequently, the vortex hamiltonian TYv obtained from is simply /C. It 
contains no "kinetic" piece and is purely "potential" . 

The symmetries of are well known and listed listed in Table 1 (see 
either [11] or [12] for a somewhat different approach). Corresponding to 
rotational invariance there is the total angular momentum 

N 

I = p,hJ2Qc.\z''\^ (10) 
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while invariance under translations gives the total linear momentum 

N N 

P = PohJ2<lc.Z°', P = PohJ2<laZ''. (11) 

a a 

P has the form of a polarization vector which justifies calling it such. It is 
also related to the total superfluid momentum density g^^ in the following 
way. Let P be a disk in the complex plane such that z'^ E T> for all a. Then 
using the identity (see Appendix) 

,a^J_ f dzAdz 

27riJ-r) z - z'^ ^ ^ 



which is independent of the size of V, we find that in rectangular coordinates 

Pj = -ejk f g'Jd^T (13) 
Jc 

where ejk is the totally anti-symmetric pseudotensor {ejk = —^kj, ei2 = 1)- 



Symmetry 


Conserved 
Quantity 


Observable 


Lagrange 
Multiplier 


Physical 
Interpretation 


Time Translation 
Invariance 


/C 


Energy 


P = l/T 


Temperature 


Space Translation 
Invariance 


P,P 


Superfluid 
Current 


E, E 


External 
Flow 


Rotational 
Invariance 


I 


Angular 
Momentum 


Q 


Total 
Vorticity 



Table 1. Table of the symmetries of the vortex lagrangian and its' conserved charges. 
The grand partition function for the ensemble is then 

Z[E,E]^ g{N,q)X{N,q)ZN[E,E] (14) 

config 

where 

Z^[E, ^] = n X '-^^^-^^ exp [-(3[1C - (EP + EP) + ^Q/]) (15) 
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{EP + EP = EjPj), is the A?"- vortex partition function and the sum is over 
all configurations of the system. g(A'", q) is the multiplicity factor due to the 
number of identical vortices in each configuration while X{N,q) represents 
the total fugacity of the system. Their exact forms are quite complicated, 
but fortunately are not needed for our purposes. Q and Ej are the lagrange 
multipliers corresponding to the conserved charges / and Pj respectively. 
Physically, Q is interpreted as the total vorticity (angular velocity) of the 
system. Since even in the absence of a net external rotational fiow in the sys- 
tem there still may be a net vorticity due to that of the vortices themselves, 
Q = ^Qa- This is the case we shall always consider. Ej is identified as the 
components of a net external non-rotational fiow or "electric" field which 
may be slowly varying in time. In complex coordinates E = {Ei +iE2)/ V^- 

First, we note that the grand partition function is a function of the 
magnitude of Ej only and not its' direction: Z[E,E] = 2[\E\] where 
\E\^ = EE = EjEj/2. This is due to the rotational invariance of JC and I. 
Second, in the absence of any external field whatsoever, {Pj)\Ej=o = for 
all values of Q; once again due to the rotational symmetry of /C and /. Such 
a state, however, is never realized experimentally and, when Q = 0, is, as we 
shall see, extremely unstable. Since in reality there is always some external 
current flow in the system, we will introduce an external flow to the system 
and look once again at the behavior of (Pj) when Ej ^ 0. The system will 
behave very differently depending upon whether or not Q vanishes and we 
will treat the two cases separately. 
3. The Case Q 0. 

Turning our attention to the N-voitex partition function, we complete 
the square by letting 

^« = ^«_i!^. (16) 
PohQ ^ ' 

Because IC is translationally invariant, 

Z[|E|] = exp(/3|E,.£;,.)z[0] (17) 



and it is now straight forward to calculate the average polarization vector 

{Pj)^mE,. (18) 
The fluctuations Pj = Pj — (Pj) about this average flow is 

{P,Pk) = (19) 
From (13), and using translational invariance in the Ej limit, 

{P,Pk)o = ne.iskm f {gtH^9n[{0))od^r (20) 

JQ 

where f2 is the total surface area of the film. As usual we define the normal 
fluid density as 

Pn^^Jjr/{r)r/mod'r. (21) 

Using (19) and (20) to calculate pn, we immediately see that the system 
will always behave as a normal fluid. There is no superfluid phase when 
Q ^ 0. Note also that for this case {Pj)\Ej=o = {Pj)o and, not surprisingly, 
rotational symmetry is restored when Ej — > 0. 
4. The Case Q = 0. 

We now sum only over those configurations of the system for which the 
net circulation vanishes. The term proportional to / is no longer present in 
the Boltzmann factor so that neither (Pj) nor {PjPk) is trivial to calculate. 
We can, however, get some notion of their behavior in the static limit by 
looking at the properties of (SH) = —{EjPj). We begin by establishing a 
very important inequality. 

Returning once again to the grand partition function, we scale the 
coordinates of the vortices by letting w" = Ez'^ so that idw^ A dw'^ = 
\E\Hdz'^ A dz"^. Then 

Z[\E\]^ V g{N,q)\{N,q)e^^(-2\N - ^D\\og\E\\zN[l] (22) 

f • \ L Pkt J / 

conftg 

where -D = X] ?a Pkt = l/cfet is the Kosterlitz — Thouless temperature. 
In obtaining this expression we have made use of the identity 

(AT \ ^ AT 

Y^Qa] =D+J2qaqp. (23) 



Then 



Pm = -{E'-^E'^yo,Z = 2{iN) 



(D) . (24) 



Since |ga| > 1 for all a, D > N and 



p{sn) < 2 1 - (iv) 




(25) 



Using this inequality and the observation that {N) > it is straightforward 
to show that {5H) has the following properties: 



1. m)\E,^o = o, 

2. (Sn) < for all P > Pkt {T < Tkt), Ej ^ 0, 

3. if {6n) = for some p > pkt {T < Tkt), Ej 0, then (N) = there, 

4. if (Sn) > for some p < pkt {T > Tkt), Ej 0, then (N) > there, 



where for completeness we have included the result obtained in Section II. 

In the absence of any vortices whatsoever, the HMW Theorem prevents 
a phase transformation from taking place and the ^He film will behave as 
a normal fluid. Since the low temperature phase is a superfluid, we will 
assume that {N) ^ when P > Pkt- From properties 2 and 3 we find that 



for all P > Pkt- A priori there is no reason why this limit will exist. If, 
however, (STi) — > — oo, as Ej 0, we would then have the unphysical result 
that an infinitesimally small external field causes an infinitely large shift in 
the average energy of the system. On physical grounds we conclude that 
{STi)o must be finite. Moreover, we will find that in order to define a normal 
fluid density we must also require that {SH) be expandable in a Taylor Series 
about Ej =0. 

We now see that due to the presence of vortices in the low temperature 
phase, {5H)\Ej=o 7^ {SH)o- Even when the applied fleld is turned off rota- 
tional symmetry remains broken so that the low temperature phase is in a 



lim (Sn) < 



(26) 



10 



state of broken symmetry. Returning for a moment to the definition of the 
superfluid velocity, we let be the phase of ip. From the definition of v^^ 
and (4), 



= -zlog 



AT 

' Z — z"- 



n 



(27) 



La=l 

where we have once again set all to zero. Now preform a global gauge 
transformation i/j — > exp [iNx)'^- Then (f) ^ (f) + Nx, or, equivalently, 
z°' — > exp {ix/2qo,)z°'. A global gauge transformation of the Bose field is the 
same as a uniform rotation of the vortices. Breaking rotational invariance is 
equivalent to breaking the U{1) gauge symmetry so that like superfluidity in 
three dimensions, superfluidity in two is also charactorized by the breaking 
of a f/(l) gauge symmetry. Furthermore, because (/C)o + {S7i)o is simply the 
average energy of the broken symmetry state, the broken symmetry state 
has less energy than the symmetric one (/C)|s^.=o and is the one favored 
energetically. 

For convenience we deflne R = —/3{5H). It is a function of the external 
held Ej so that by expanding i? in a Taylor Series, 



R — Rq + 



2! dEkdEk 



EjEj + H.O.T. (28) 





There are no terms in the expansion which containing odd powers of Ej since 
^ is a function of |£'| only. The coefficient of each term in the expansion, 
although independent of Ej , is still a function of /5 and the vortex chemical 
potentials. Moreover, because of (26) the expansion starts with a constant 
term below the Kosterlitz — Thouless temperature. We can now formally 
solve for the grand partition function in terms of the applied field 

Z[|^|] = (y^)''"exp(A)Z (29) 

where Z is a function of P and the chemical potentials only, while 



A=- 



1 1 d^R 



2 2! dEkdEk 
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EjEj + H.O.T. (30) 



Note that ^ — > when £'j — > 0. We are now in the position to look at the 
behavior of the correlation functions. 
The average momentum is simply 

/?(P,) = i?o^ + ||- (31) 

which diverges as 1/\E\ when Ej below Tkt- Notice, however, that 
P{S7i)o — —Rq is still finite. The fluctuations about this inflnite average 
current in the Ej limit is 



^ \ J k/o ^^^^ ^^^^ y jk ^^^^ J 



6jk. (32) 





2 dEidEi 

Taking the trace of (32) we find that because we are working in two dimen- 
sions 

(33) 



"^'^ 2m(5VLdEidEi 
We see that pn is finite as long as R is expandable in a Taylor Series about 
Ej = 0, justifying our anzatz that {STC)o exists. Then, using (13), we 
obtain the following expression for the superfluid momentum-momentum 
correlation function 

P'ejiekm Jj9i^{r-)rJ{0))od'r= Jim^ (^6,^ - 2^) + Pmpr,S,k 

(34) 

where Ug = Ro/Vl. 

When rig 7^ we flnd that the off-diagonal terms in the superfluid 
momentum-momentum correlation functions are inflnite when Ej ^0, sig- 
nifying off-diagonal long range order in the static limit. Since the high 
temperature phase is a normal fluid, must vanish above the transition 
temperature Tg although it is greater than zero below it. functions as 
an order parameter for the system. It is straightforward to show that this 
order parameter must also be independent of (3 below T^. Calculating the 
total average energy of the system and using (29), we obtain the following 
consistency equation 

R IdRo OA IdZ 
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Now let Ej — > 0. Whether or not (/C)o is finite in this hmit is determined 
not by how the hmit is taken, but rather by the expression for /C itself. Z 
is independent of Ej and, from (30), A ^ as Ej ^ 0. Because R is finite 
in the Ej — > limit, we find that Rq must be independent of /3 everywhere 
except, perhaps, at the transition temperature where the grand partition 
function itself is no longer analytic. The order parameter is therefore a step 
function in T with the discontinuity occurring at Tc- It is, however, still a 
function of the chemical potentials of the vortices. 

From property 2 we conclude that > T^t- Prom our inequality 
(25), we note that although it may be well behaved everywhere else, there is 
singularity in (A^) when T = Tkt- If Tc > Tkt, the system will have a singular 
behavior before the transition temperature, which would be unphysical. On 
physical grounds, then, we conclude that Tc = Tkt- 

There is not much more that we can say about the general ensemble. 
We now turn our attention to the Kosterlitz — Thouless ensemble. Of course 
all the results we have obtained so far also holds for the Kosterlitz — Thouless 
ensemble with the addition of a few new features. Using the form of the 
grand partition function given in (22) to calculate j3'^{PjPk)o directly, we 
find that 



The inequality is now exact since D = N ior the Kosterlitz — Thouless en- 
semble. Moreover, 



where N = N — (N)^^ is the fluctuation in the average number of vortices. 
It vanishes in the static limit for all values of T except, perhaps, at the 
transition temperature T^t- Because there are now only vortices with charge 
|g| = 1, there is only the one fugacity A and, by using (37), we conclude that 




(36) 




(37) 



d{N)!^ 



kt 



A 



dX 



= 0. 



(38) 
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(A'")q* is independent of the chemical potential and from (36) we find that 
for the Kosterlitz — Thouless ensemble the order parameter n^* is a constant 
independent of all thermodynamic variables. 

We now let Ej vary slowly with time and look at how the average 
number of vortices in the system changes. From (33) and (37) we find that 
for small Ej, 

d{N)^' _ 1 pj'n ds 

dt Tkt -T ks dt ^ ' 

where 

E = \mEjEj (40) 
is the amount of energy the external source is depositing into the system. 

■ kt 

Notice that (A^) is coupled only to the normal fluid density and not to the 
superfluid density. Above the transition temperature the external source 
will decrease the average number of vortices while below the transition tem- 
perature it will tend to increase it. 

Returning to (36), we find that the average density of vortices for the 
Kosterlitz — Thouless ensemble p^* = {N)q^/Q, in the static limit is 

kt rp 

Pv=^jr^(^iTkt-T) (41) 

where 9{x) is the step function. Consequently, there are no vortices in the 
high temperature phase, a result which runs contrary to the standard model 
of the transition and with which we shall attempt to reconcile at the end of 
this paper. Rather, the number of vortices is unbounded at T^t and slowly 
decreases as T is lowered below T^t. Because of the attractive force between 
oppositely charged vortices, this decrease is due to the annihilation of a 
charge q = 1 vortex with a. q = —1 vortex. This does not mean, however, 
that the fluid is quiescent at high temperatures. In our previous work [10] 
we have demonstrated that the excitation spectrum of the ^He fllm consists 
of two components: one a quasi-particle component corresponding to the 
phonon gas, and the other a psuedo-particle, vortex gas component. In 
the constant density limit, the vortex gas component completely decouples 
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from the phonon gas so that the interaction between the quantized point 
vortices and the quasi-particle spectrum may be neglected to lowest order. 
This is the regime that we are working in. Consequently, there will in fact 
be quasi-particle excitations in the fluid in the high temperature phase, and 
these excitations may even consist of rotational flows. Our result states that 
these flows may not form ideal, point vortices with integer valued circulation. 

Let us return for a minute to the General Kosterlitz — Thouless ensem- 
ble. Because the energy required to create higher charged vortices is much 
greater than the energy need to create vortices with unit charge, < p^* 
where is the average density of vortices in the general ensemble. This 
observation, combined with the inequality (25), gives an upper bound on 
the vortex density in the general ensemble 

n T 

P-^^jr^O{nt-T) (42) 

as well as < n^* (see flgure 1). The average number of vortices in the 
general ensemble may decrease faster with temperature than the number 
of vortices in the Kosterlitz — Thouless ensemble. How much faster is still 
an open question, although we can obtain a formal expression for it. From 
(24), we find 

P. = -y^T^^(7^fct-T)+7^ (43) 



where 



^=^T^^^^^°"^^^"^ ^^^^ 



is the approximate difference in the rates at which the average number of 
vortices decrease. Note also that because 7?. — > as T — > 0, {D)q — {N)q ~ 
{T/TktY where p > 0. 
5. The Behavior at T = Tu. 

We now turn our attention to the behavior of the system at the tran- 
sition temperature itself. While before we fixed (5 ^ Pkt and considered the 
behavior of the correlation functions as Ej — > 0, we now fix Ej at some 
small but non- vanishing value and let /? — > /?^^ (T — > T^). Our results in 
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this section holds for both the general ensemble as well as the Kosterlitz — 
Thouless ensemble. At the transition temperature the system must be scale 
invariant. Kosterlitz and Thouless have shown that the relevant scale for 
the ensemble is the cutoff length a and to emphasize this fact we shall, in 
this section, denote the grand partition function as Z[a, Let us now 
scale a ^ ^a. Then 



ZNi^a, \E\\ [[ — -— ^ — exp I (jckt QaQpU 

(45) 

Letting — z"- we find that Z[^a, \E\] = Z[a,$,\E\] so that scaling a is 
equivalent to scaling the external flow fleld. Thus, the system must also be 
invariant under the scaling of Ej at the transition temperature. 

Returning to (29) and using the definition of the normal fiuid density, 
we write 

^ n Q 1 
logZ[a, \E\] = logZ + ^ log (EkEk) + -m(3p^nEkEk + H.O.T. (46) 

Leaving the logarithmic term alone for now, we see that the normal fiuid 
density is the coefficient of a term which has naive scaling dimension 2. 
As such, we conclude that ^ as T — > T^. Following Minnhagen and 
Warren (see Note Added), we define the superfiuid density as ps = Po~Pni 
and find that 

Ps_ _ Po_ _ 2mk i 

where is Boltzmann's constant and we have used ksTki = e^t- This is the 
result originally obtained by Nelson and Kosterlitz [3] (our ps is the number 
density and not the mass density used by Nelson and Kosterlitz, however). 
Moreover, we find that all the higher order terms in the expansion of R, 
which are proportional to the higher order superfiuid momentum correlation 
functions, must also vanish as T ^ T^. 

Ug is the cofficient of the logarithmic term in (46) which has an anoma- 
lous scaling dimension. It is unclear exactly what, if any, restrictions scale 
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lim ^ = ^ = f^ (47) 



invariance will impose upon it. If, however, we take the scaling invariance 
to also include the logarithm, we find that ns{T = Tkt) = which sub- 
stantiates our earlier physical arguments for the vanishing of Ug above the 
transition temperature in Section IV. 
6. Conclusion. 

From the above results we are lead to a somewhat different interpreta- 
tion of the Kosterlitz — Thouless superfluid phase transition. Approaching 
the phase transition from above we see that there are no vortices whatsoever 
above the transition temperature. From the HMW Theorem the ^He film 
has no choice but to behave as a normal fluid. As T is lowered below the 
Kosterlitz — Thouless temperature, however, the vacuum state becomes un- 
stable and any external perturbation of the system causes a gas of vortices 
to be created. Due to angular momentum conservation, this gas must be 
neutral and consequently the ^He film behaves as a superfluid. The HMW 
Theorem is circumvented due to the presence of vortices in the fluid. As 
the temperature is lowered further, the average number of vortices in the 
system decreases due to pairwise annihilation of oppositely charged vortices 
until there are no vortices left at T = 0. Approaching the transition tem- 
perature from below, the increasing temperature creates vortices pairwise 
in the fluid to preserve neutrality. The average separation between vortices 
now increases with temperature, thereby decreasing the likelihood of pair 
annihilation. This seperation, however, is bounded above by the finite size 
of the system while the average number of vortices will continue to grow 
without bound. At the transition temperature itself a massive annihilation 
of the vortices will occur so that no vortices are left above Tj^f 

At flrst glance this seems to be different from the standard model of 
the phase transition which has free vortices above the transition tempera- 
ture and vortices bound in dipole pairs below it. We note, however, that in 
the Kosterlitz — Thouless theory two vortices are considered to be bounded 
when they come within a distance a, the cut-off length, of each other. From 
our point of view they have effectively annihilated each other. What Koster- 
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litz and Thouless interpreted as vortex binding and unbinding may also be 
viewed as vortex annihilation and creation. Our p^* is, in the language of 
the standard model, the density of free vortices in the liquid and long range 
order is destroyed by these vortices when the number of vortices in the sys- 
tem becomes infinite in a finite size system. Aside from the requirement that 
there are no vortices in the high temperature phase, there is little difference 
between our theory and the Kosterlitz — Thouless theory. 

Although we have allowed Ej to vary with time, ours is not a complete 
theory of the dynamical phase transition because the vortices have not been 
coupled to the underlying fluid. Ej is treated as an external field and has 
not been related to any of the elementary excitations in the ^He film. The 
reader is referred to either [13] or [14] for the complete description of the 
dynamical phase transition. 

As we have seen, the magnitude of the external field |£'| plays an anal- 
ogous role to the chemical potential for the system. It enters, however, into 
the grand partition much in the same way as the external magnetic field 
does for the two dimensional Ising Model; as a lagrange parameter. Identi- 
fying as the chemical potential, equation (29) shows that when Ug > 
there is a zero of the partition function at = 0. This is reminiscent of the 
Yang — Lee theory of phase transitions in which a phase transition occurs 
when the partition function develops a zero as the complexified fugacity z 
pinches the real axis at z = 1 [15], [16]. The crucial difference, of course, is 
that Ej is a vector, and not a scalar like the chemical potential. 

Note Added 

There are two different definitions of the superfluid density currently in 
use in the literature. The one used by Nelsen and Kosterlitz is 





while the one given by Minnhagen and Warren is 
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where we have used our definition of g^^ verses Minnhagen and Warren's 
(they difi^er by a factor of ph/m). Minnhagen and Warren has shown that 
the two definitions are equivalent near T^f. Our definition of ps agrees 
with Minnhagen and Warren's once the difference of the factor of two in 
the definitions of the vortex hamiltonian (8) is taken into account. All three 
definitions of the superfiuid density give the same result at Tkt, namely (47). 
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APPENDIX 

Due to the importance of the relationship between the polarization 
vector and the superfluid velocity, we shall establish the identity (12). We 
begin with a very brief discussion of differential forms. The reader is referred 
to [17] for a complete description. 

The object A in the integration measure idzAdz defines a wedge product. 
Its' usefulness will be made clear later on. Let / be any function of the 
coordinates Xj. The exterior derivative df of / is defined as 

df = ^dxj {Al)z 

with the rule dxj A dxk = —dxk A dxj. In evaluating the integral in 
(12), we find it much more convenient to use polar coordinates and we 
let z = rexp{i6). Then from (Al), dz = {dr + irdO) ex];){i6) so that 
dz Adz = —2irdr A d9 which is obtained without the use of the Jacobian of 
the transformation. Then 

2m Jt) z-z'^ ttz'^ Jo Jq \ z^ J 



where R is the radius of the disk V. We have broken the integral up into 
two pieces so as to make use of the expansion 

^ CX) 

l3^ = E-^ (^3) 
which holds as long as \z\ < 1. Because 

r e'^'d9 = 6n,o, (AA) 



we find that 

1 r dzAdz 2 
2m Jjy z- z'^ 

and we are done. 



— / rdr = z"" {A5) 
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FIGURE CAPTIONS 



Figure 1. Graph of the density of vortices p^* verses T for the Kosterhtz — 
Thouless ensemble. The shaded region in the graph represents the allowed 
values that the density of vortices in the general ensemble may have. No- 
tice that there are no vortices whatsoever above the transition temperature 
and the singularity at T = T^t- 



